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ABSTRACT 

This paper describes a very general approach to the calculation of the Zeeman split- 
ting effect produced by an external magnetic field on the rotational levels of diatomic 
molecules. The method is valid for arbitrary values of the total electronic spin and 
of the magnetic field strength -that is, it holds for molecular electronic states of any 
multiplicity and for both the Zeeman and incomplete Paschen-Back regimes. It is based 
on an efficient numerical diagonalization of the effective Zeeman Hamiltonian, which 
can incorporate easily all the contributions one may eventually be interested in, such 
as the hyperfine interaction of the external magnetic field with the spin motions of the 
nuclei. The reliability of the method is demonstrated by comparing our results with 
previous ones obtained via formulae valid only for doublet states. We also present re- 
sults for molecular transitions arising between non-doublet electronic states, illustrating 
that their Zeeman patterns show signatures produced by the Paschen-Back effect. 

Subject headings: magnetic fields — polarization — molecular processes — methods: 
numerical 



Introduction 



Most polarized radiation diagnostics of astrophysical magnetic fields have been carried out 
via the theoretical interpretation of the observed polarization signatures in atomic spectral lines 
(e.g., the reviews by Bagnulo 2003; Mathys 2002; Stenflo 2002). However, over the last few years 
we have witnessed an increasing interest in molecular spectropolarimetry as a tool for empirical 
investigations on solar and stellar magnetism, concerning both the molecular Zeeman effect (e.g., 
the recent overviews by Asensio Ramos & Trujillo Bueno 2003; Berdyugina et al. 2003 and Landi 
Degl'Innocenti 2003a; see also Uitenbroek et al. 2004) and the Hanle effect in molecular lines (e.g., 
Landi Degl'Innocenti 2003b; Trujillo Bueno 2003a). 
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The aim of the present paper is to describe in some detail our approach to the Zeeman and 
Paschen-Back effects in diatomic molecular lines, which we have been applying over the last few 
years while contributing to the development of the field of radiative transfer in molecular lines 
(Asensio Ramos & Trujillo Bueno 2003, 2005; Asensio Ramos et al. 2003, 2004a, 2004b, 2005). 
Our strategy is based on a a very efficient numerical diagonalization of the effective Zeeman Hamil- 
tonian, whose general expression has been derived by a number of authors working in the field 
of molecular physics (e.g., Brown & Carrington 2003). Interestingly, in physics laboratory exper- 
iments, where the applied magnetic field is known beforehand, the observation of the splittings 
provide a measurement of the magnetic dipole moment of the considered molecule in the particular 
spin-rotational level involved. Given that the magnetic moment depends on the electronic struc- 
ture of the molecule, it is obvious that its measurement provides information about the molecular 
structure. In astrophysics we have the inverse problem, the magnetic field being the unknown 
quantity. To obtain information about cosmic magnetic fields, therefore, we have to make use of 
our most precise knowledge on molecular structure in order to infer the magnetic field vector via 
the theoretical modeling of polarization signals in molecular lines. To this end, we need: 

• to obtain the molecular number densities at each point within the astrophysical plasma model 
under consideration, 

• to calculate the splittings of the molecular energy levels with the strengths of the individual 
Zeeman components, and 

• to solve the radiative transfer problem for the emergent Stokes parameters. 

Here we focus on the issue of calculating the splittings and strengths of molecular transitions, 
including the most general situation of the incomplete Paschen-Back effect for electronic states 
of arbitrary multiplicity. For information on the techniques we use for calculating the molecular 

number densities and for solving the Stokes- vector radiative transfer problem in (magnetized) stellar 
atmospheres we refer the reader to the papers by Asensio Ramos et al. (2003) and Trujillo Bueno 
(2003b), respectively. 

It is of historical interest to mention that the Zeeman effect in diatomic molecules was consid- 
ered shortly after the development of the quantum theory. Kronig (1928) investigated the molecular 
Zeeman effect in Hund's (a) and (b) cases for the angular momentum coupling between electronic 
and rotational motion. Only one year after. Hill (1929) investigated the Zeeman effect for doublet 
states of diatomic molecules in intermediate states between Hund's cases (a) and (b). The status 
of the theory was reviewed by Crawford (1934), emphasizing that, at that moment, the Zeeman 
effect was understood in pure Hund's cases (a) and (b) and in intermediate cases between the two. 
At that time, the intensity of the Zeeman transitions had not been investigated in detail because 
the calculations were rather involved when using the basis functions of Hund's case (b). Fifty 
years after the paper of Kronig, Schadee (1978) re-investigated the Zeeman effect for doublet states 
of diatomic lines in the intermediate case between Hund's cases (a) and (b), but using the basis 
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functions of Hund's case (a). This greatly facilitated the calculation of the intensity of the Zeeman 
transitions. The issue of the Zeeman effect in lines of diatomic molecules was later considered by 
Tiling (1981) who studied in great detail Schadcc's (1978) theory and applied it for understanding 
the broad-band circular polarization observed by Harvey (1973) in near-IR lines of CN. 

All the above-mentioned developments were based on a number of approximations for the 
description of the molecular motions for the zero field case. They take into account the rotational 
energy (but neglecting centrifugal distortions) and the strongest angular momenta couplings. For 
states without electronic orbital angular momentum (S states), only the spin-rotation coupling 
was included. For states with non-zero electronic orbital angular momentum (11, A, . . . ), only the 
spin-orbit coupling was taken into account. 

The formulae developed by Schadee (1978) are only applicable to doublet states of diatomic 
molecules in the Zeeman or Paschen-Back regimes. Recently, Berdyugina & Solanki (2002) have 
extended Schadee's (1978) formulation to allow for the calculation of the effect of a magnetic 
field on states with arbitrary spin, but limited to the Zeeman regime. Their strategy consisted 
in numerically diagonalizing a simplified Hamiltonian corresponding to the zero-field case and 
accounting for the Zeeman Hamiltonian as a first order perturbation. This approach, which neglects 
the non-diagonal matrix elements (A J 7^ 0) of the Zeeman Hamiltonian, is only valid in the linear 
Zeeman regime. Therefore, the maximum reliable value of the magnetic field strength is established 
by the transition to the Paschen-Back regime. 

In this paper we present a very general approach which allows us to calculate the effect of 
a magnetic field on the rotational levels of diatomic molecules in electronic states with arbitrary 
multiplicity. The method is valid in both the Zeeman and the Paschen-Back regimes. It is based on 
the numerical diagonalization of the effective molecular Hamiltonian, which describes the molecular 
motion using the basis functions of Hund's case (a). Therefore, the inclusion of any additional 
contribution to the effective Hamiltonian of the diatomic molecule is straightforward, provided the 
matrix elements of the Hamiltonian are known in the basis functions of Hund's case (a). This makes 
it possible to investigate effects like hypcrfine structure (HFS) without much additional effort. Such 
refinements are not only of interest in physics laboratory experiments, but also when it comes to 
interpreting correctly the linear polarization signals that anisotropic radiation pumping processes 
induce in spectral lines (e.g., Landi Degl'Innocenti &; Landolfi 2004). The illustrative examples 
shown in this paper neglect HFS effects, but we plan to consider this interesting molecular HFS 
problem in a future investigation. 

The outline of this paper is the following. Section 2 and the various appendices describe 
our numerical diagonalization approach in some detail, including a summary of some results from 
angular momentum theory with the aim of facilitating a better understanding of the paper. In order 
to verify the reliability of our numerical results, in Section 3 we compare them with those that can 
be obtained via Schadee's (1978) theory. We will show also some applications to non-doublet states. 
Finally, Section 4 summarizes our main conclusions with an outlook to future research. 
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2. Angular momentum theory for diatomic molecules 

The theory of angular momentum provides a robust theoretical framework for describing the 
complex structure of diatomic molecules. In this section we summarize how the angular momentum 
theory and Racah's tensor algebra allow the description of the quantum mechanical behavior of 
any diatomic molecule under the presence of an arbitrary magnetic field. Although the theory on 
which our work is based on can be found explained in detail in specialized monographs on molecular 
spectroscopy (e.g., Judd 1975, Brown & Carrington 2003), we have decided to present here a brief 
summary in order to facilitate the understanding of this paper and to avoid any possible notational 
confusion. 

2.1. Introduction 

As is well-known, the Born-Oppenheimer approximation leads to an effective separation of 
the energies due to the electronic, vibrational and rotational motions in molecules (e.g., Herzberg 
1950). This approximation is supported by the fact that the mass of the nuclei is several orders of 
magnitude larger than the mass of the electrons. This implies that, when nuclei move, the electrons 
adapt rapidly to the new nuclear configuration. As a result, the electrons feel some kind of effective 
potential, which depends only on the positions of the nuclei and on the electronic configuration 
of the particular electronic state. The Born-Oppenheimer approximation inmediately leads to the 
possibility of separating the eigenfunctions associated with each motion. 

For convenience, we have selected Hund's case (a) eigenfunctions because they lead to several 
simplified expressions (see Appendix A for details). In particular, the eigenvalues of the angular 
momentum operators are very easy to obtain in this basis set, which greatly simplifies the calculation 
of the matrix elements of the total Hamiltonian. The obvious consequence of selecting a basis set, 
which is a truly good basis only in some limiting conditions, is that the total Hamiltonian is usually 
non-diagonal and it has to be diagonalized to obtain the energies of the molecular levels. In the 
usual notation, the Hund's case (a) eigenfunctions can be written as (Judd 1975): 

\aASJ:;v;njM) = \aAST.)\v)\njM), (1) 

where we have explicitly indicated that the total eigcnfunction is a product of the eigenfunctions 
associated with the electronic, vibrational and rotational motions (which is only strictly valid under 
the Born-Oppenheimer approximation). 

We point out that \aAST,) is the electronic eigcnfunction. The symbol a represents a collection 
of quantum numbers which are used to label the electronic configuration of the molecule, while S 
is the total spin. The symbols A and S are the projections of the total orbital electronic angular 
momentum (L) and of the total spin (S) on the internuclear axis, respectively. Both are good 
quantum numbers in Hund's case (a), while S is not a good quantum number in Hund's case (b) 
(see Appendix A for more details). The vibrational eigenfunction is represented by \v). Because 
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we are interested in lines of diatomic molecules, this eigenfunction can be correctly described with 
only one quantum number v. Finally, \QJM) represents the rotational eigenfunction. This function 
depends on the total angular momentum J, its projection M on the quantization axis (typically 
chosen as the axis along the magnetic field vector) and the number Q = |A| + S. It is interesting to 
note that $7 is a good quantum number only in Hund's (a) case, since S is only well defined in this 
case. Note also that, due to the presence of the absolute value, in the non-rotating molecule the 
energy of the spin-orbit levels with J7 > 1/2 are doubly degenerate. When the molecule is rotating, 
this degeneracy is broken (see, e.g. Herzberg 1950). A brief summary of the good quantum 
numbers in Hund's cases (a) and (b) can be found in Table 1. Additionally, vector diagrams of the 
two coupling cases are shown in Figure 11, which clearly show the strength of each coupling. 

Due to the difference in mass between the electrons and the nuclei of the diatomic molecule, 
it is appropriate to refer the motion of the electrons to a frame F' fixed to the nuclei whose origin 
is at the center of mass of the molecule, rather than to an external laboratory frame F (with 
its origin also at the center of mass of the molecule). Fig. 1 illustrates the relative position of 
both frames. The coordinates of an electron in the laboratory frame F are {x,y,z), while the 
coordinates on the frame fixed to the molecule F' are (^, rj, Given the arbitrariness in the choice 
of the frames, it is convenient to consider the ( axis of the F' frame parallel to the intcrnuclcar axis 
and the z axis of the F frame along the quantization axis. This quantization axis is usually chosen 
along the direction of the magnetic field vector when it is present. Obviously, in the absence of a 
magnetic field, it can be arbitrarily chosen in any direction. Since both reference frames share the 
same origin, coordinates in one frame can be transformed to the other frame by means of a simple 
rotation. This rotation can be parameterized in terms of the three Eulcr angles (a, /?, 7). Once the 
coordinates of an electron are known in one of the frames, a standard geometrical rotation between 
both frames using the Euler angles can be used to obtain the coordinates in the other frame. As 
usual, once the ( axis of the F' frame is fixed, there is an additional freedom in the choice of the 
direction of the ^ and 77 axis, which depends on our choice for the Euler angle 7. Without loss of 
generality, we choose 7 = 7r/2. 

The normalized rotational eigenfunctions \CIJM) in Hund's case (a) can be obtained by solving 
the Schrodingcr equation for a symmetric top. The solution depends explicitly on the Euler angles, 
and can be written as (Judd 1975) 



where P][^j^(a, /3, 7) is the rotation matrix (e.g., Edmonds 1960). The representation of the rota- 
tional eigenfunctions in terms of the rotation matrices is very convenient because it facilitates the 
subsequent calculations. A fundamental property is that of conjugation for the rotation matrices, 
which can be expressed as 




(2) 



Vi,a{a,l3,jr = {-l) 



'-M-n(«,A7)- 



(3) 



Other important properties of the rotation matrices which are used in this section are summarized 
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in Appendix B. Note that the conjugation property ahows us to verify that the rotational eigen- 
functions are orthonormal. In fact, one can apply Eq. (3) together with Eq. (B2) to verify that 
{n'J'M'\njM) = Sii'nSj'jSM'M- 

The rotation matrices are of great help for transforming tensorial operators from one frame 
to another, provided that the Euler angles between both frames are known. Given that we are 
dealing with rotations, it is advantageous to consider the spherical components of a tensor instead 
of working with the cartesian ones (see, e.g., Edmonds 1960 for the definition of the spherical 
components of a tensor). Let us consider in some detail the case of a vector, since it will be helpful 
for calculating the matrix elements of the Hamiltonian. Assume that r is a vector in whose 
spherical components in the laboratory frame F are given by fg, where q = 0, ±1. The spherical 
components of the same vector in the frame fixed to the molecule F' are denoted by r^. The 
relationship between both components can be written as (e.g., Edmonds 1960): 

r<i= Yl rgVl,g{-a,-(3,-^), (4) 

q'=0,±l 

where we have used the fact that the rotation needed to carry the frame fixed to the molecule to 
the laboratory frame is expressed by the inverse rotation (—a, — /3, —7). Using the properties of the 
rotation matrices (see, e.g., Edmonds 1960), the previous equation can also be written as: 

E rgVl^'{(^,l3,^r. (5) 

q'=0,±l 

The final part of this introduction to the mathematical tools needed to calculate the matrix 
elements of the Hamiltonian is related to the Wigner-Eckart (WE) theorem (e.g., Edmonds 1960). 
This theorem facilitates the calculation of the matrix element of any tensor operator because it 
takes full advantage of any symmetry that may be inherent to the problem under consideration. 
In other words, it isolates those parts of a problem that are essentially geometric in character from 
those which depend explicitly on the physics of the problem. For the q component of a given tensor 
of rank k, the WE theorem reads: 

(JM|rf|j'M') = (-l)'^-^(^_^ J ^^yjWT^'^J'). (6) 

The quantity (J||r('^) || J') is called the reduced matrix element, a number which depends on the 
physics of the selected problem. An example is the reduced matrix element that involves an angular 
momentum J between its own eigenfunctions: 

(J||J|| J) = VJ(J + 1)(2J + 1) (7) 

2.2. Molecular Hamiltonian 

In this section we present the total effective Hamiltonian which we use to describe the molec- 
ular motion including the effect of a magnetic field. Consider a diatomic molecule in a given 
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electronic, vibrational and rotational state and in the presence of an external magnetic field. It is 
possible to describe the molecular motion by writing the full Hamiltonian which takes into account 
all the Coulomb forces among electrons and nuclei and the Lorcntz force due to the presence of 
a magnetic field. However, it is more appropriate in terms of both economy and feasibility to 
consider the effective Hamiltonian approach (Brown et al. 1979, Brown & Carrington 2003). This 
effective Hamiltonian has the same eigenvalues as the full Hamiltonian, but it operates only within 
the rotational, spin and hyperfine levels of a given vibrational level belonging to a given electronic 
state. The non-diagonal matrix elements (coupling the vibrational level of interest with any other 
vibrational level of any other electronic state) arc absorbed in the effective Hamiltonian and repre- 
sented by some parameters which can be obtained by comparison with laboratory or astrophysical 
observations. Following Brown &; Carrington (2003) (see also Brown et al. 1979) we write down 
the following very general expression for the effective Hamiltonian: 

HeS = Hso + Hss + -f^rot + ^cd + ^sr + -f^LD 

+ -f^cdLD + -f^Ms + -f^cdhfs + Hz- (8) 

The term Hso represents the coupling between the total spin of the molecule S and the orbital an- 
gular momentum L. H^^ represents the coupling between the spins of the electrons in the molecule. 
The rotational energy of the molecule is accounted for by the term Hj-ot, while i?cd represents the 
contribution to the rotational energy of the centrifugal distortion as a consequence of the rotation 
of the molecule. The interaction between the total molecular spin and the rotational angular mo- 
mentum is described by iJgr- The additional degeneracy present for the two possible projections of 
L on the internuclear axis, namely ±|A|, is broken by the inclusion of the terms Hj^d and -f^cdLD 
of the effective Hamiltonian. These terms represent the A-doubling interactions (interaction with 
higher energy electronic states) and the appropriate centrifugal distortion corrections, respectively. 
In case one is interested in taking into account the hyperfine structure, it suffices to include the 
terms i?hfs and i?cdhfs) which are the hyperfine interaction and the corresponding centrifugal distor- 
tion of the hyperfine interaction as the molecule rotates, respectively. Such hyperfine contributions 
will not be considered in this paper, but they will be the subject of a future investigation. Finally, 
Hz represents the interaction between the molecule and an external magnetic field B. 

Each one of the terms of the Hamiltonian can be written with the aid of the different angular 
momentum operators that can be defined in the molecule. As previously mentioned, it is useful 
to write all these operators using the spherical tensor notation because it leads to considerable 
simplifications in the notation and it allows the use of the powerful tools of the angular momentum 
theory. The same notation we have used for vectors can be straightforwardly extended to the 
spherical components of tensorial operators. On the one hand, will refer to the q spherical 
component of the tensor A*^ of rank k in the molecular fixed frame F'. On the other hand, the 
spherical component of the same tensor in the laboratory frame F will be indicated by A^. This 
extension to tensors of arbitrary rank k is only needed when we include in the Hamiltonian tensorial 
operators (see, e.g., Brown &; Carrington 2003). In the effective Hamiltonian used in this paper 
we consider only vectorial operators with rank up to A; = 1. However, it is important to emphasize 
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that this formalism allows us to include any contribution to the effective Hamiltonian written in 
terms of tensorial operators, which constitutes one of its big advantages. 

The explicit form of each of the terms included in the effective Hamiltonian has been obtained 
previously by molecular physicists (see Brown & Carrington 2003, and references therein). The 
effective Hamiltonian approach is flexible enough to allow writing the explicit form of some of the 
terms of the Hamiltonian in different ways. It is advantageous to write them in a way that simplifies 
the calculation of the matrix elements in the chosen basis set. The spin orbit coupling term reads 

Hso = ALlS^ + ^An [N^LIS^ + L^S^f^^] , (9) 

where A and Ad are the spin-orbit coupling constant and the centrifugal correction to A, re- 
spectively. The first part includes the coupling between the total spin and the orbital angular 
momentum. Note that this contribution using Hund's case (a) basis set is just proportional to 
the product of A and S since these are the Lq and Sq components of L and S, respectively. The 
second term includes the effect of the centrifugal distortion in the spin-orbit coupling. The symbol 
stands for the scalar product of N by itself, i.e., = N • N. Following Brown & Carrington 
(2003), wc use as the square of the rotational operator instead of (see Appendix A). The 
latter approach is preferred by some authors because this operator is proportional to the rotational 
kinetic energy. However, using R as the rotational operator involves the angular momentum L 
which has non-diagonal elements between different electronic levels. In the effective Hamiltonian 
approach, the effect of these non-diagonal terms are absorbed into the rotational and coupling 
constants. This way, we end up with a Hamiltonian which operates only within the (rotational, 
spin and hyperfine) levels of a given vibrational level of an individual electronic state. 

Focusing now on the spin-spin interaction, several functional forms have been suggested in 
the literature (e.g.. Brown & Carrington 2003 for a summary). It has been shown that they are 
equivalent and we have decided to choose the following functional form which greatly simplifies the 
calculation of the matrix elements in Hund's case (a) basis set: 

i^ss = |a [3{Sl)' - S2] , (10) 

where A is the spin-spin coupling constant. Obviously, this contribution turns out to be important 
only when more than one electron is uncoupled (S > 1/2). 

Turning our attention to the rotational part of the effective Hamiltonian, the pure rotation 
Hamiltonian can be written as 

Hrot = BN^, (11) 

where B is the usual rotational constant For the evaluation of the matrix elements in the basis set 
of Hund's case (a), it is better to represent the rotational angular momentum N in terms of J and 
S since these operators have well-defined eigenvalues in this basis set. This transformation can be 
easily carried out because N = J — S (see Appendix A). Concerning the inclusion of centrifugal 
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distortion terms, it is sufficient to arrive up to the third order contribution even for laboratory 
experiments, and the exphcit form of the Hamiltonian turns out to be: 

H^d = -D{N^)'^ + F(N2)3, (12) 

where D and H are the quartic and sextic distortion constants. Finally, the coupling of the spin 
and the rotation can be included in the effective Hamiltonian as the scalar product between the 
rotational angular momentum and the spin: 

H,r = 7(J - S) • S, (13) 

where 7 is the spin-rotation coupling constant (note that this constant is different from the 7 Euler 
angle). In this paper we will not include more terms in the effective Hamiltonian for the zero- field 
case. 

An additional constraint has to be included when dealing with diatomic molecules that greatly 
simplifies the evaluation of the matrix elements of the effective Hamiltonian. This constraint is 
related to the fact that the molecule does not rotate around the internuclear axis and so (Jq — 
Lq — Sq) = (where the component of the rotational angular momentum operator is along the 
internuclear axis due to the selection of the frame F'). 

When the molecule is under the action of a magnetic field B, the magnetic sub levels pertaining 
to each energy level of total angular momentum J have slightly different energies, essentially due 
to the precession of the total angular momentum around the quantization z-axis, which we have 
chosen along the magnetic field vector. It is possible to write a Zeeman effective Hamiltonian for a 
diatomic molecule that takes into account the coupling between the magnetic field B and the total 
spin S and orbital angular momentum L and between the magnetic field and the rotation of the 
molecule N = J — S: 

Hz = gsfJ-o^ ■ S + gLlJ-oB ■ L 

- 5,/xoB ■ (J - S), (14) 

where /Lto is the Bohr magneton, B is the magnetic field vector, gs is the electron spin ^-factor, 
is the electron orbital gi-factor and gr is the rotational gi-factor (including the nuclear and electronic 
contribution g^. = g^^^^^^^ — gf^^). If the spin of the nuclei is non-zero, a contribution of the coupling 
between the nuclear spins and the magnetic field can be included with: 

Hr = Y.9hmB-I\ (15) 
1 

where the nuclei are labeled by i. Because we neglect hyperfine structure, we will also neglect this 
term. 

Although small differences exist, it is usual to put 55 = 2 and = 1- The exact value of gs is 
slightly larger than 2 (its value can be obtained from Quantum Electrodynamical considerations) 
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and in fact the values of the ^-factors can be slightly different due to the inclusion of additional 
perturbations from other electronic states. On the contrary, due to the difference of mass between 
the electron and the proton, the rotational (j(-factor (which is of the order of /i7v//iO) with ^jsj the 
nuclear magneton) is ~3 orders of magnitude smaller than the electronic ^-factors. Therefore, the 
contribution of the last term of the Zeeman Hamiltonian is only relevant when the molecule has 
no spin and no orbital angular momentum. In this case, there is no contribution to the Zeeman 
effect coming from the electronic motion because the electronic cloud is essentially spherical and 
the molecule is in a state. This is the case of the fundamental electronic states of CO and SiO, 
for which the sensitivity to the presence of a magnetic field is very small. Ideally, the value of 
the ^-factors should be determined by confronting the theoretical description of the molecule with 
spectroscopic observations. The rotational ^-factor includes a pure nuclear contribution and another 
one coming from the coupling between the rotation of the electronic cloud and the magnetic field 
(Judd 1975). The nuclear contribution to the rotational (/-factor for a diatomic molecule formed 
by two nuclei A and B can be estimated by the following formula (Judd 1975): 

nuclear _ MiV {Za^II + ZuM^) 

^'^ MO MaMb {Ma + Mb) ' 

where Ma and Mb are the masses of each nuclei in atomic mass units and Za and Zb are the 
atomic number of each nuclei. It is important to take into account that the importance of the 

rotation of the electronic cloud can be decisive. For example, the calculated value for the nuclear 
contribution to the OH rotational (/-factor is g^'^'^^^'^^ ~ 5.25 x 10~^ while the measured value is 
gr = -6.33 X 10"^ (Brown et al. 1978). 

A simplification of the functional form of the Zeeman hamiltonian arises when the quantiza- 
tion axis z lies along the magnetic field direction. In this case, the scalar products between the 
magnetic field vector and the angular momentum operators in Eq. (14) are transformed into their 
corresponding projections along this axis multiplied by the modulus of the magnetic field vector. 
These projections are the spherical components of the operators along the z axis of the laboratory 
frame: 

Hz = gslJ-oBSl + gLUoBLl - griioB{J^ - Sq). (17) 

2.3. Matrix Elements 

The energy levels of a diatomic molecule in the presence of an arbitrary magnetic field can be 
obtained via the numerical diagonalization of the total effective Hamiltonian matrix. In contrast to 
previous approaches, ours does not make use of any perturbational calculation since we diagonalize 
the full Hamiltonian including the Zeeman contribution. The first step is to write the effective 
Hamiltonian matrix in the chosen Hund's case (a) basis set. The advantage of this basis set is 
that the majority of the molecular fixed frame components of the angular momentum operators are 
diagonal in this basis, which largely simplifies the evaluation of the Hamiltonian matrix elements. 



(16) 
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Of course, for problems where Hund's case (a) is a good representation of the angular momentum 
coupling of the particular rotational level of the electronic level under consideration, the Hamil- 
tonian will be highly diagonal. In general, the Hamiltonian will have a non-diagonal contribution 
that manifest deviations from the Hund's case (a) coupling. 

One of the easiest ways to calculate the matrix elements of the effective Hamiltonian is to sub- 
stitute the rotational eigenf unctions by their explicit expressions in terms of the rotation matrices 
given by Eq. (2). The complex conjugate of the rotational eigenfunctions (jr2M| are obtained by 
taking the complex conjugate of the rotation matrix and transforming them with the help of the 
conjugation property given by Eq. (3). Since many of the matrix elements of the angular momen- 
tum operators are known in the molecular fixed frame, it is convenient to evaluate them in this 
frame. In case the components of a tensor in the laboratory frame appear in the Hamiltonian (for 
example in the Zeeman Hamiltonian), a transformation to the molecular fixed frame is performed 
using the standard transformation rules for the tensorial components given by Eq. (5). Finally, the 
calculation of the matrix elements of a spherical component of an arbitrary operator Ag is reduced 
to an integration over the Euler angles. The integrands of the {J'^' M'\Ag\JQM) integrals consist 
in the product of two or three rotation matrices (two from the rotational eigenfunction and an extra 
one from the transformation to the molecular fixed frame). These integrals are easily calculated 
using the Weyl's theorem given by Eq. (Bl) and the property given by Eq. (B2) in Appendix B. 
In other cases, it is advantageous to directly use Wigner-Eckart's theorem given by Eq. (6). For 
facilitating the calculation and the reproducibility of our results, we give the analytical expressions 
for the matrix elements of each term of the Hamiltonian in Appendix C. 



2.4. Diagonalization 

From an inspection of the explicit formulae for the matrix elements of the effective Hamiltonian 
in Appendix C, we can see that many of the terms included for the description of the diatomic 
molecule are diagonal in several quantum numbers. For instance, the Hamiltonian associated to 
the spin-orbit coupling is completely diagonal in Hund's case (a) basis set. The molecular rotation 
contribution is diagonal in A, S, J and M, but not in^ S and CI, while the Zeeman interaction is 
diagonal only in A, S and M. Consequently, in the non-zero field case, the total angular momentum 
J is not a good quantum number. Therefore, only M remains as a good quantum number and the 
effect of the magnetic field has to be described in the Paschen-Back regime. Strictly speaking, J is 
not a good quantum number for any arbitrary value of the magnetic field strength. Obviously, for 
sufficiently weak magnetic fields, J is an approximately good quantum number and we can safely 
treat the molecule in the Zeeman regime. 



^These are good quantum numbers in the case of a non-rotating molecule, but they are not good quantum numbers 
when rotation is taken into account. 
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When selecting the basis set, we have to make sure that the chosen set of lASE; QJM) functions 
has to be complete enough to capture all the non-diagonal terms in the Hamiltonian matrix. In 
our case, the Hund's case (a) basis set was chosen using the following rules: 

• We include eigenfunctions with the two possible values of the projection of the orbital angular 
momentum on the internuclear axis A, i.e., ibA. This is fundamental to appropriately take into 

account the effect of A-doubling in which there is a break of the degeneracy between energy 
levels with the same value of [A|. Since in this paper we do not include any A-doubling 
interaction, the eigenvalues corresponding to both values of A arc degenerate. However, if A- 
doubling needs to be explicitly accounted for, it is only a matter of including the appropriate 
effective Hamiltonian and calculating its matrix elements in Hund's case (a) basis set. 

• We include eigenfunctions with the two possible values of the projection of the spin on the 
internuclear axis S, i.e., ibS. 

• When calculating the energy of the magnetic sublevels of a given level with total quantum 
number J, we also include the effect of the non-diagonal terms between the J level and the 
levels with J — 1 and J+1. This is only necessary when a magnetic field is present because the 
contribution of the Zeeman splitting is the only term included in the effective Hamiltoninan 
which is non-diagonal in the quantum number J. In some cases (for instance, for heavy 
molecules), one may be interested in including non-diagonal terms with A J = ±2. With the 
present approach, this is only a matter of augmenting the basis set by taking into account 
J — 2, J— 1, J, J+1 and J + 2 for a given level J. 

• For each value of J, we include all the possible values of M from — J to J. It is important 
to note that the Hamiltonian is always diagonal in M unless the hyperfine structure is taken 
into account. 

The total number of eigenfunctions included in the basis set and, consequently, the size of 
the effective Hamiltonian matrix, is obtained by summing all the possible values of the quantum 
numbers: 2 possible values of A (1 in case A = 0), 25 + 1 possible values of S and 2J + 1 values of 
M for each included value of J. Therefore, the size of the basis set can be expressed in terms of 
the spin of the electronic state and the J value of the level of interest: 



For instance, let us assume we choose a ^H electronic state. The eigenfunctions of the form 
ASTi-jQJM) used for the calculation of the energies of a rotational level J in this state are 



1, 1/2, 1/2; 3/2, J',M), |1, 1/2, -1/2; 1/2, J',M), |-1, 1/2, 1/2; -1/2, J',M) and |-1, 1/2, -1/2; -3/2, J',M) 



(we have separated each number with commas for clarity reasons). Since we include the coupling 
for the values of J between J — 1 and J + 1, we have to build such a set of eigenfunctions for each 




3(25 + 1)(2J + 1) if A = 0,J>1 
6(25+ 1)(2J+ 1) if A7^0,J>1. 



(18) 
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value of J' = J — 1, J, J + 1. Additionally, we take into account the values of M for each allowed 
value of J' that are in the range — J' to J'. Summing up all the cigenfunctions obtained from the 
previous combinations, we end up with 24 J + 12 eigenf unctions, a number which can be quite high 
even for small values of J. 

The Hamiltonian matrix built with these eigenfunctions is a J\f x J\f matrix. Due to the 
orthogonality properties of the eigenfunctions of Hund's case (a), many of the matrix elements 
are zero. As a result, the Hamiltonian matrix is very sparse (the number of non-zero elements 
is usually much less than the number of zero elements). In order to numerically diagonalize this 
matrix, we can make use of any of the available numerical procedures (see, e.g.. Press et al. 1986), 
although algorithms specifically built for sparse matrices should be used in this case. As stated 
in Eq. (18), the size of the matrix increases linearly with the value of J and the calculation of 
the eigenvalues and eigenvectors of a big matrix represents a very hard numerical work even for 
intermediate values of J. However, we have seen that the effective Hamiltonian matrix is diagonal 
in the subspaces spanned by eigenfunctions with different values of M (or Mp in case we include 
the hyperfine structure). In this case, if we reorganize the matrix by ordering the basis set by 
their M quantum number, we end up with a block-diagonal matrix. Each block belongs to the 
space spanned by the eigenfunctions with a given value of M. Therefore, we transform the problem 
of diagonalizing a J\f x J\f matrix to the diagonalization of 2 J + 1 submatrices of size 6(25 + 1). 
This is shown in Figure 2, which corresponds to an example with J = 2. The standard algorithms 
applied to calculate the eigenvalues of a general matrix scale as 0{J\f^), so that this reorganization 
translates into a huge decrease in the computing time (of the order of (2 J + 1)^). 

With the numerical diagonalization procedure, we calculate the eigenvalues and eigenvectors 
of the effective Hamiltonian matrix. The eigenvalues are associated to the energies of the magnetic 
sublevels M of a given rotational level of the electronic state of interest. The energy shift produced 
by the presence of a magnetic field can be obtained by the difference between these energies and the 
zero-field energies, obtained by diagonalization of the effective Hamiltonian neglecting the Zeeman 
Hamiltonian Hz- The eigenvectors can be used then for calculating the expectation value of any 
operator. In our case, we are interested in calculating the expectation value of the dipolar moment 
operator, that is related to the relative strength of the transitions. The eigenfunction associated 
with any level can be written as the following linear combination of the eigenfunctions of Hund's 
case (a): 



where the sum is extended over the values of A, S, Q and J included in the basis set. Since the 
matrix elements of the dipolar moment operator are known in Hund's case (a) (sec, e.g., Schadee 
1978), it is possible to use the previous linear combination to calculate it in the new basis. 




(19) 



Asnj 
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3. Illustrative Results 

This section is devoted to showing some comparisons between the results obtained via our 
numerical diagonalization of the Hamiltonian and with the previous approach to the molecular 
Zeeman effect based on Shadee's (1978) theory. Although the formalism of Schadee (1978) allows 
one to calculate the Zeeman splittings and strengths of the Zeeman components of transitions 
among doublet states (5 = 1/2), it has been recently extended to transitions among states with 
arbitrary spins but for the linear Zeeman regime (Bcrdyugina & Solanki 2002). Our approach is 
suitable also for calculating the Zeeman splittings and strengths of the components for the case 
in which a transition to the Paschen-Back regime occurs for electronic states with arbitrary spins. 
The Hamiltonian diagonalization approach is the most general one, since it permits us to include 
easily any additional coupling among the angular momenta of the molecule and is valid for arbitrary 
strengths of the external field. Moreover, it can be applied to any electronic state with an arbitrary 
value of the spin. 

Since one of the various applications we are carrying out is the synthesis of Stokes profiles 
induced by the molecular Zeeman effect in strongly magnetized regions of the solar atmosphere, wc 
show first the Zeeman patterns for several lines of different diatomic species which are observed in 
the solar atmosphere. As is well known, a Zeeman pattern diagram shows the position and relative 
strength of each of the a (AM = ±1) and tt (AM = 0) components which arise from transitions 
among the magnetic sublevels of the upper and lower rotational levels of a given spectral line. The 
components in the upper part of the diagram are the a components, with the AM = — 1 transitions 
indicated as vertical lines going upwards and the AM = 1 transitions indicated as vertical lines 
going downwards. The relative strength of each component is proportional to the length of the 
vertical line. Additionally, the tt components are in the lower part of the diagram. We have 
selected several molecular lines from MgH, OH, CN, C2, which show interesting Zeeman patterns 
(see, e.g., Schadee 1978; Illing 1981; Berdyugina k, Solanki 2002; Asensio Ramos & Trujillo Bucno 
2003; Asensio Ramos et al. 2005). In addition, we show also the Zeeman patterns for the CCS 
radical, which is a molecule of interest in the research field of star formation regions. Since this is 
a linear molecule, it can be described with the same formalism as if it were a diatomic molecule. 

3.1. Doublet states 

3.1.1. MgH 

Our first example is the P2(5.5) line of the A^Yl - (0,0) electronic band of MgH. The 

ensuing MgH lines, which are located in the visible range around 5150 A, have recently become of 
interest because, apart from presenting an interesting magnetic sensitivity, they show conspicuous 
scattering polarization signals when observed close to the solar limb (Stenflo & Keller 1996, 1997). 

Figure 3 shows the Zeeman patterns obtained for a magnetic field strength of 1000 G. The 
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left panel concerns the results obtained using the theory developed by Schadee (1978) while the 
right panel shows the results obtained using our numerical diagonalization of the Hamiltonian. The 
molecular constants have been obtained from Hubcr &; Herzbcrg (2003). For consistency with the 
assumptions of Schadee (1978) and in order to properly compare the results, we have included 
the very same terms in the effective Hamiltonian. In this case, the total Hamiltonian of the A^H 
electronic state has only contributions from the spin-orbit coupling, while that of the lower X^T,'^ 
level has only contributions from the spin-rotation coupling. Figure 3 shows that both Zeeman 
patterns are indistinguishable, a proof that our numerical diagonalization technique is working 
properly. As previously indicated by Schadee (1978) (see also Berdyugina & Solanki 2002), such 
MgH lines present strong Paschen-Back effects. The reason is that the lower electronic state has a 
very small spin-rotation constant and the energy gap between two consecutive energy levels is very 
small. Therefore, the presence of a weak magnetic field leads to a Zeeman splitting comparable 
to this energy gap, which produces interferences between the magnetic substatcs. The presence of 
Paschen-Back effects produce a deformation of the Zeeman pattern with respect to the symmetric 
situation of the linear Zeeman regime. The ultimate cause of such interferences is that the Zeeman 
Hamiltonian Hz is not diagonal in J, so that the energy levels do not have a J value associated 
with them. When a magnetic field is present each energy level has contributions from levels having 
different values of J in the zero-field case. 



3.1.2. OH 

The next two examples shown in Figures 4 and 5 are for the Pi(10.5) and P2(9.5) near-IR 
lines of OH, which belong to the vibro-rotational (2,0) band of the fundamental electronic state 
X^H. Here the transition to the Paschen-Back regime occurs for magnetic strengths slightly below 
10^ G (Berdyugina & Solanki 2002). Therefore, for magnetic fields of ~1000 G we are safely in 
the Zeeman regime. Similarly to the MgH case, the left panels show the results obtained with the 
theory of Schadee (1978), while the right panels show the results that we have obtained when using 
the numerical diagonalization of the Hamiltonian including only those terms which were considered 
by Schadee (1978). The molecular constants have been also obtained from Huber & Herzberg 
(2003). Note the extremely good agreement between both calculations. Since these OH lines are 
in an intermediate coupling scheme between Hund's case (a) and (b) (Berdyugina & Solanki 2002), 
this calculation constitutes another proof that our numerical diagonalization technique is working 
properly. Although we have written the Hamiltonian matrix using the basis functions of Hund's 
case (a), we can correctly recover the behavior for levels which are not correctly described under 
this case. In this situation, the Hamiltonian is non-diagonal and the final eigenfunctions result 
from appropriate combinations of the Hund's case (a) eigenfunctions. 

The circular polarization profiles produced by these OH lines in sunspot umbrae were first 
observed by Harvey (1985) and modeled by Berdyugina & Solanki (2001) taking into account that 
the effective Lande factor of the two lines have opposite signs (Riiedi et al. 1995). The fact that 
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both lines have effective Lande factors of opposite sign can be noted in Figs. 4 and 5. To this 
end, we recall the definition of the effective Lande factor as the wavelength shift from line center of 
the center of gravity of the AM = —1 components in Lorentz units (i.e., in units of = XqUl/c, 
where Aq is the central wavelength, c is the speed of light and vr, is the Larmor frequency of the 
magnetic field). The Zeeman patterns show that the displacement from line center of the AM = — 1 
components have opposite signs for both lines. 

With the polarization sensitivity we can obtain nowadays with the available solar polarimeters, 
the effective Hamiltonian we have used appears to be sufficient for the description of the OH 
transitions. However, we must emphasize that additional effects, like for instance the A-doubling, 
can be easily included in case a more realistic description turns out to be needed for explaining the 
emergent Stokes profiles. 

3.1.3. CN lines in the IR 

Another interesting example is the calculation of the Zeeman patterns for a molecular transition 
with high J-values, since it permits to test the numerical diagonalization approach for a case in 
which the effective Hamiltonian matrix is very large. One example is the Qi(46.5) near-IR line of 
CN observed by Asensio Ramos et al. (2005). This line belongs to the A^H — X^S+ electronic 
transition of CN and it is located in the same spectral region as the two OH lines discussed in the 
previous example. The rotational constants were again taken from Huber &: Herzberg (2003). The 
magnetic field at which the transition to the Paschcn-Back regime occurs for the levels of the upper 
A^n electronic state is ~560 kG (e.g., Berdyugina & Solanki 2002). Therefore, we can safely treat 
them in the Zeeman regime for the typical stellar magnetic fields. However, the levels of the A^H 
electronic state have to be described in an intermediate coupling scheme between Hund's cases (a) 
and (b). On the other hand, the levels of the lower electronic state X^E"*" can be correctly described 
under the Hund's case (b) coupling when the magnetic field is weak enough so that we are in the 
linear Zeeman regime. However, the transition to the Paschen-Back regime in the lower electronic 
state occurs for rather weak fields (~77 G for the lowest J levels). For this reason, these lines are 
always in the Paschen-Back regime under the typical magnetic fields of sunspots. This example 
represents then a complicated problem in which neither the Zeeman regime nor any of the limiting 
Hund coupling cases can be applied. 

Asensio Ramos et al. (2005) have pointed out that the observed polarization signal from 
this and other similar lines present an anomalous behavior, which is due to the peculiar form of 
the Zeeman patterns. In Fig. 6 we show the Zeeman patterns for this line calculated via the 
Hamiltonian diagonalization technique. We have verified that these results are similar to those 
obtained using the formulation of Schadee (1978). Figure 6 shows the patterns for weak fields 
(50 G), intermediate fields (500 and 2500 G) and very strong fields (30000 G). The first one 
shows that, even for weak fields, the CN lines cannot be correctly described under the Zeeman 
regime and that interferences between close J levels are of importance. At the typical kG field 
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strengths of sunspots, the AM = ±1 components have almost the same center of gravity while 
that of the AM = components has a displacement with respect to the AM = ±1 components. 
This peculiarity, produced by the Paschen-Back regime, leads to the explanation of the anomalous 
emergent Stokes profiles in the umbral spectrum. 

3.1.4. CN lines in the UV 

Another case of interest concerning transitions between doublet states is the B^T,^ — X^S"*" 
ultraviolet transition of CN. This band has recently gained some attention due to the strong linear 
polarization signals observed close to the solar limb (Stenfio 2003, Gandorfer 2003). Theoretical 
investigations of the scattering polarization and the Hanlc effect in these molecular lines have been 
carried out by Asensio Ramos &: Trujillo Bucno (2003, 2005) in order to explain the ladder structure 
of the observed Q/I signals. Since the magnetic field in the observed quiet regions is assumed to be 
relatively weak and the scattering polarization calculations are very demanding, these theoretical 
investigations were performed by using the Lande factors obtained under Hund's case (b) coupling. 

Here we show that the Paschen-Back effect in this electronic system is of importance for fields 
of the order of 1000 G (and even for fields below 100 G). The rotational constants were taken 
from Huber Sz Herzberg (2003). We have plotted only the results obtained via the numerical 
diagonaUzation of the total Hamiltonian, but they are similar to those obtained using the formulae 
of Schadcc (1978). Figure 7 shows the Zeeman patterns for the i?i(70.5) line for a field strength 
of 1000 G. This result is representative of the Paschen-Back effect for lines between high J levels. 
The iii(70.5) line is one of the spectral lines presenting strong scattering polarization signals in the 
spectral region between 3771 A and 3775 A. It is interesting to note that Paschen-Back effects are 
important for fields as low as 1000 G for transitions between high- and low- J levels. Therefore, for 
the typical kG fields of sunspot umbrae the lines of this system should present anomalous Stokes 
profiles. The observable effects of these perturbed Zeeman patterns on the emergent Stokes profiles 
should be similar to those observed in the lines of the ^4^11 — system, and whose Zeeman 

patterns are shown in Figure 6. 

3.2. Non-doublet states 

3.2.1. C2 

The first case we have considered is the (0,0) transition of the electronic band d^U — a^U of 
C2; i.e., one of the so-called Swan band. The investigation of this band in the solar atmosphere has 
gained some interest recently due to the linear polarization signals produced by scattering processes 
(see, e.g., Gandorfer 's 2000 atlas of the linearly-polarized solar limb spectrum) and their sensitivity 
to the presence of a weak magnetic field via the Hanle effect (e.g., Trujillo Bueno et al. 2004). In 
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this paper we are however interested in calculating the Zeeman patterns of some lines of this band 

for a magnetic field of 3000 G. We have chosen the Pi (7) and P2(6) lines as representative of the 
results which can be obtained also for other lines of the Swan system. The results we have obtained 
via our numerical diagonalization of the total Hamiltonian arc shown in Figure 8. 

The critical magnetic field at which one starts to detect signatures of the transition to the 
Paschen-Back regime for this electronic system is larger than 70000 G (Berdyugina & Solanki 
2002). Since we are interested in much weaker fields, the magnetic sensitivity of the lines of this 
transition can be safely modeled in the Zeeman regime. This can also be inferred from the very 
symmetric appearance of the Zeeman patterns shown in Figure 8. Furthermore, the deviations 
from Hund's case (b) are not very large and we have verified that the obtained Zeeman patterns 
have some similarities with those obtained in Hund's case (b). An important conclusion regarding 
the effective Lande factor can be obtained by a quick look to the Zeeman patterns. Note that the 
center of gravity of the AM = — 1 transition is very close to zero in the P2 (6) line while it is much 
larger for the Pi (7) line. This shows that the Pi lines are much more magnetically sensitive than 
the P2 lines. The same is valid for the P3 lines. 

It is very difficult to find molecules having electronic states with S ^ 1/2 that present strong 
Paschen-Back effects and which can be observed in solar and/or stellar spectra in the visible and 
infrared regions. The target molecular transitions are those in which the energy difference between 
consecutive levels is small enough so that the Zeeman splitting for typical stellar magnetic fields 
produce significant interferences between levels with consecutive values of J. Since the spin-orbit 
coupling is usually much larger than the spin-rotation one (see Huber & Herzberg 2003 for the 
molecular constants), a candidate molecule should ideally have an electronic state with A = so 
that the small spin-rotation coupling gives strong interferences between closely lying energy levels. 
A few molecules observed in the solar and other stellar atmospheres present such states. This is 
the case of CO with its a'^T,~^ electronic state, which is radiatively linked with the fundamental 
X^T,~^ level by a forbidden transition around 1800 A. 

Another case of interest is C2, with the Ballik- Ramsay system. This system results from the 
electronic transition — a^n„ situated in the infrared (Ballik Sz Ramsay 1963), with the center 
of the band around 1.8 /xm. As an example, we show in Figure 9 the Zeeman patterns for the Ri(16) 
and Pi (16) lines. Both lines present strong Paschen-Back effects. Curiously, the deformed Zeeman 
patterns resemble those patterns we calculated for the CN lines of the A^H — X^'E'^ system. In fact, 
we expect that such C2 lines present Stokes profiles that resemble those observed in sunspot umbrae 
for the previously considered CN lines. Thus, they would present very weak circular polarization 
signals and stronger and antisymmetric linear polarization signals. Such C2 transitions can be easily 
observed in carbon stars (Goorvitch 1990) and they have been also detected in comets (Johnson & 
Larson 1983). In this respect, it is of interest to mention that Johnson & Larson (1983) give a list 
of possible observable bands in comets. The (0,0) band is around 5632.7 cm"-*^, the (1-0) band is 
around 7080.8 cm~^ and the (2-0) band is around 8506.5 cm~^. 



3.2.2. SO and CCS 



Turning now our attention to the millimeter-wave region, it is worthwhile to mention that it 
is possible to observe pure rotational transitions in the lower rotational levels of molecules with 
ground electronic states having A = and S ^ 1/2. Some of these molecules are present in regions 
of star formation and the investigation of the Zeeman effect in these molecules is very important for 
gaining information about the magnetic field. This is the case of SO and the CCS radical, which are 
used as tracers of dense cores in star formation regions. Since they present no hyperfine structure, 
and due to the observed narrow and intense emission peaks, they have been used to measure the 
velocity structure in dense cores (Langer et al. 1995). Due to the two unpaired electrons present 
in both molecules, the Zeeman splittings are important, and this fact has been used for searching 
for magnetic fields in the densest parts of star formation regions. Shinnaga & Yamamoto (2000) 
have calculated the Zeeman effect in the fundamental electronic state of both SO and CCS in 
order to give accurate Lande factors for several rotational lines. Figure 10 demonstrates that we 
are able to reproduce their Fig. 2 for the Ri(3) transition of CCS. Of course, we have used the same 
rotational and coupling constants in order to be able to compare properly both results. The results 
of Shinnaga & Yamamoto (2000) were obtained by diagonalization of the effective Hamiltonian 
written using Hund's case (b) eigenfunctions. In this case, it is more appropriate because, for 
weak rotation, the angular momenta coupling is closer to Hund's case (b). Remarkably, using our 
approach based on Hund's case (a) eigenfunctions, we obtain the same results, thus demonstrating 
again that our approach is very robust. Interestingly, we find that we need to use a magnetic field 
strength of 1 mG instead of 100 /xG in order to reproduce their results. Apparently, the magnetic 
field indicated in the labels of Fig. 2 in Shinnaga & Yamamoto (2000) appears to be inconsistent 
with the results given in the text, that are otherwise correct. Finally, it is important to note that, 
since the energy separation of rotational levels decreases as the weight of the molecule increases, it 
might be possible that the AJ = ±2 matrix-elements are non-negligible for a correct modeling of 
the Zeeman splitting in some molecular species. 



4. Conclusions 



The numerical approach presented in this paper allows us to calculate the effect of a magnetic 
field on the energy levels of diatomic molecules for arbitrary values of the total electronic spin and of 
the magnetic field strength (i.e., it is valid for states of any multiplicity and for both the Zeeman and 
incomplete Paschen-Back regimes). The ensuing computer program we have developed gives the 
splittings of the molecular energy levels and the strengths of the individual Zeeman components 
(i.e., the Zeeman patterns of molecular transitions), which is the basic information we use in 
our radiative transfer code for modeling the emergent Stokes parameters from magnetized stellar 
atmospheres. It is based on an efficient numerical diagonalization of the effective Hamiltonian, 
which gives the eigenvalues and eigenvectors for each magnetic sublevel, from where we obtain the 
energies of the molecular levels and the expectation values of the dipole moment operator between 
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the pair of levels producing the molecular transition under consideration. 

Our numerical diagonalization approach of the effective Hamiltonian can treat cases in which 
quantum interferences between very close in energy levels lead to the transition to the Paschen- 
Back regime. We take into accoTint the effect of non-diagonal terms between levels with A J = ±1 
in a self-consistent way, without making use of any perturbativc calculation. Our strategy gener- 
alizes previous ones which only allowed the investigation of doublet states for arbitrary values of 
the magnetic field (Schadee 1978) or states with arbitrary spin but in the linear Zeeman regime 
(Berdyugina & Solanki 2002). Additionally, our computer program allows a straightforward inclu- 
sion of any additional term into the effective Hamilonian, so that the description of the molecular 
motion can be as refined as needed. 

We have performed several comparisons between our results and those obtained via previous 
formulations based on Schadee's (1978) theory, demonstrating that both results are indistinguish- 
able when the very same terms of the effective Hamiltonian are included. We have also performed 
calculations for several molecular transitions arising from non-doublet states (e.g., states with 
S = 1). Since the theory developed by Schadee (1978) is limited to doublet states, we have per- 
formed comparisons in the Zeeman regime with the results given by pure Hund's coupling cases. 
Furthermore, we have also calculated Zeeman patterns for diatomic lines in triplet states in which 
clear Paschen-Back effects are at work. 

It has been clearly shown in physics laboratory experiments that a very detailed description of 
the molecular motion is needed to investigate correctly the magnetic properties of diatomic molec- 
ular lines (e.g.. Brown & Carrington 2003, and references therein). However, in astrophysics the 
investigation and application of the Zeeman and Hanle effects in molecular lines is still at an early 
stage of development. Fortunately, a new generation of polarimeters and telescopes should allow us 
to obtain high-quality spectropolarimetric observations with unprecedented spectral resolution and 
polarimetric sensitivity. We believe that the numerical diagonalization approach described in this 
paper constitutes the ideal framework to investigate in detail the magnetic properties of molecular 
lines in different astrophysical environments, since it permits a rather straightforward inclusion of 
all the desired terms in the effective Hamiltonian. Among the various investigations we have in 
mind for the near future, we would like to mention the influence of the hyperfine structure and/or 
the A-doubling on both the Zeeman and Hanle effects in the spectral lines of diatomic molecules. 

We thank Egidio Landi Degl'Innocenti and the anonymous referee for their careful reading of 
our paper and for suggesting some useful improvements. This research has been funded by the 
European Commission through the Solar Magnetism Network (contract HPRN-CT-2002-00313) 
and by the Spanish Ministerio de Educacion y Ciencia through project AYA2004-05792. 
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Fig. 1. — The two frames used to describe the motion of a diatomic molecule. The F frame with axes 
(x, y, z) is fixed at the center of mass of the molecule, with the z axis along the quantization axis. 
The F' frame with axes t], (") rotates with the molecule, with the ( axis along the internuclear 
axis. The transformation of any tensorial quantity between both frames can be easily carried out 
with the aid of Eqs. (4) and (5). 
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Fig. 2. — Rotational part of the effective Hamiltonian matrix for a level with J = 2. A similar 
matrix has to be built for any combination of A and E. We indicate with a "x" the elements which 
can, in principle, be different from zero. The rest of matrix elements are zero due to the orthonor- 
mality of the eigenfunctions. Since we are indicating only the rotational part of the Hamiltonian, 
the eigenfunctions are of the form \JM), with $7 prescribed by the values of A and S. The left 
panel shows the Hamiltonian matrix when the basis set is sorted by the value of J, that results in 
a highly non-diagonal matrix. The right panel shows the same Hamiltonian matrix when the basis 
set is ordered by the value of M. Since the effective Hamiltonian is diagonal in the M quantum 
number, the resulting matrix is block-diagonal, thus allowing a block diagonalization in each of the 
subspaces spanned by the eigenfunctions with the same value of M. 
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Fig. 3. — Zeeman patterns for a magnetic field strength of 1000 G for the P2(5.5) fine of the MgH 
electronic transition A^H — The left panel shows the results applying the theory developed 

by Schadee (1978) while the right panel shows those obtained from the numerical diagonalization 
of the Hamiltonian. This low- J transition of MgH is in the transition to the Paschen-Back regime 
for a field of 1000 G, thus the Zeeman patterns are highly perturbed due to the presence of the 
non-diagonal AJ = ±1 terms in the total effective Hamiltonian. 
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Fig. 4. — Zeeman patterns for a magnetic field strength of 1000 G for the P2(9.5) line between 
vibrational levels of the fundamental X^Il state of OH. The left panel shows the results applying the 
theory developed by Schadee (1978) while the right panel shows those obtained from the numerical 
diagonalization of the Hamiltonian. The magnetic field at which there are interactions between 
adjacent rotational levels in this electronic state is much higher than 1000 G, so that this line is in 
the Zeeman regime. As a consequence, the Zeeman patterns are almost unperturbed. However, the 
transition has to be described using an intermediate coupling between Hund's cases (a) and (b). 
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OH P, (10.5) - 1000 G 
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Fig. 5. — Same as Figure 4 but for the Pi (10.5) line between vibrational levels of the fundamental 
X^n state of OH. 
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CN Qi (46.5) - 3QQ00 G 




Fig. 6.— Zeeman patterns for 50, 500, 2500 and 30000 G of the Qi(46.5) CN line. Note that the 
(7 and TT components tend to be sjnnmetric for low fields and get deformed due to the transition to 

the Paschcn-Back regime. At very high fields, the symmetry is again recovered. This is an example 
of a transition in which the AJ = ±1 matrix elements in the effective Hamiltonian are non-zero 
even for fields as low as 50 G. 
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Fig. 7. — Zeeman patterns for a magnetic field strength of 1000 G for the Ri(70.5) hne of the 
ultraviolet electronic transition S^S"^ — X^S"^ of CN. Note that the Paschen-Back effect is present 
even for high J values. 
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Fig. 8. — Zeeman patterns for a magnetic field strength of 3000 G for the Pi (7) and P2(6) lines of 
the electronic transition d^Il — a^U. of C2. In this case, since 5 = 1, the theory of Schadee (1978) 
cannot be applied and these results have been obtained with the numerical diagonalization of the 
Hamiltonian. 
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Fig. 9. — Zeeman patterns for a magnetic field strength of 1000 G for the Ri(16) and Pi (16) lines 
of the electronic transition — a^Iiu of C2 (the Ballik-Ramsay system). This is an example of 
a transition in which one of the levels has S ^ 1/2 while having A = 0. Since the energy separation 
of the multiplet levels is small, the AJ = ±1 matrix elements of the effective Hamiltonian are 

of importance because the splitting for moderate magnetic fields are of the order of the energy 
separation. The Zeeman patterns clearly show this perturbation in the transition to the Paschen- 
Back regime. 
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CCS Ri(3) - 1 mG 



-600 -400 -200 200 400 600 

Splitting [Hz] 

Fig. 10. — Zccman patterns for a magnetic field strength of 1 niG for tlie Ri(3) pure rotational line 
of CCS in the fundamental electronic state. This line has to be described in an intermediate 
coupling case between Hund's case (a) and (b). The comparison with the results of Shinnaga 
&; Yamamoto (2000) indicate that our approach based on a Hund's case (a) basis set works for 
transitions which are almost perfectly described under a pure Hund's case (b) coupling. 
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Case (a) 


v>_yasc I u I 


Good quantum numbers 


A, S, J, E (fi = A + E) 


A, S, J, N 


Degeneracy (non-rotating) 


\ 


f 2 ifA^O 

1 if A = 


( 2(25" +1) ifA/0 
\ 25 + 1 if A = 


Degeneracy (rotating) 


\ 


f 2 ifAy^O 
1 ifA = 


(2 ifA^O 
\ 1 ifA = 


Values of J 


|n|,|n| + i,|n| + 2,... 


\N-S\,\N-S\ + 1,...,N + S 


Values of N 


Not defined 


A,A + 1,... 


Conditions 


AA > BJ 


AA < BJ 



Table 1: Brief description of Hand's cases properties. 

A. Hund's cases 

This Appendix describes briefly the most common angular momentum coupling cases present 
in diatomic molecules. All angular momentum vectors of the molecule (electronic orbital L, spin 
angular momentum S and rotational angular momentum R) form a resultant total angular mo- 
mentum which is designated by J. Hund's coupling cases are distinguished by the strength of the 
coupling among all angular momenta present in the molecule. 

In Hund's case (a), illustrated in the left panel of Figure 11, the orbital angular momentum 
is strongly coupled to the internuclear axis by electrostatic forces. The spin angular momentum is 
in turn strongly coupled to the orbital angular momentum through spin-orbit coupling. Addition- 
ally, the rotational motion is weakly interacting with the spin and orbital motions. In this case, 
the projection of the total electronic angular momentum ft is well defined and composed of the 
projection of the orbital angular momentum L on the internuclear axis (A) plus the projection of 
the spin angular momentum S on the internuclear axis (S). The rotational angular momentum R 
thus couples with the angular momentum along the internuclear axis fl to form the total angular 
momentum J. The condition for the suitability of Hund's case (a) is that the spin-orbit coupling 
AAT, has to be much larger than BJ{J + 1), with B the rotational constant and A the spin-orbit 
coupling constant. 

In Hund's case (b), illustrated in the right panel of Figure 11, the spin angular momentum S is 
very weakly coupled to the internuclear axis. In this case, the orbital angular momentum L is still 
strongly coupled to the internuclear axis with projection A. Because of the weak coupling between 
the spin and the internuclear axis, it is not possible to define ^l. The vector A couples then to 
the rotational angular momentum R to form the resultant N (the total angular momentum apart 
from spin). Finally, the angular momentum N couples with the spin S to form the total angular 
momentum J. The condition under which Hund's case (b) applies is that the spin-orbit coupling 
must be much smaller than BJ{J + 1). 

Table 1 gives a summary of the good quantum numbers in both Hund's cases (a) and (b), 
together with information on the degeneracy of the levels in the rotating and non-rotating molecule. 
We also give the possible values of J and N (if defined) and the conditions under which each Hund's 
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case applies. 



B. Useful properties of the rotation matrices 

This appendix gives some important and useful properties of the rotation matrices which can 
be found also in angular momentum books (see, e.g., Edmonds 1960, Brink & Satchler 1968, Judd 
1975). Two properties arc used for deriving the matrix elements of the Hamiltonian. The first one 
is the Weyl's theorem, which states that the integral of the product of three rotational matrices 
over the Euler angles can be easily calculated as follows: 

/ da dPsmpvi^^^^{a,(3,j)V-'^^^^{a,(3,j)vi}^^^{a,/3,^) = 

Jq Jo Jo 

= Stt^ I "^^ "^^ "^^11 "^^ "^^ I fBl) 

\ Mi M2 Ms J yni n2 J ' ^ ' 

where we have introduced the 3-j symbols according to their standard definition (see, e.g., Edmonds 
1960). The other property of interest is 



/•27r /"TT /'27r 

i da d^ d/3sin/?P4^^(a,/3,7)*P]^^^^(a,/3,7) 
Jo Jo 



2J + 1 

where 6ab is the Kronecker's delta, which is 1 for a = 6 and zero otherwise 



'^JiJ2'^MiM2fcin2) (B2) 



C. Matrix elements of the effective Hamiltonian 

This appendix gives the explicit form of the matrix elements of the effective Hamiltonian in 
Hund's case (a) basis set. One of the purposes is to show the explicit form of the matrix elements 
in order to recognize whether or not a contribution to the total Hamiltonian is diagonal in Hund's 
case (a) basis set. As it is mentioned in this paper, these expressions have been obtained by 
applying the Wigner-Eckart theorem to evaluate the tcnsorial form of the effective Hamiltonian 
given in Section 2.2, together with an extensive use of Racah's algebra and the properties outlined 
in Appendix B (see Brown &; Carrington 2003 for details). If any of the terms is not diagonal in 
one of the quantum numbers of Hund's case (a) basis functions, we indicate it by using primed and 
non-primed quantities for the bra and the ket of the matrix element, respectively. For instance, the 
rotational Hamiltonian given by Eq. (C3) is diagonal in A, S, J and M, but not in E and ft. 



C.l. Spin-orbit coupling Hamiltonian 



{ASYinjM\Hso\AS^njM) = ^AE + AdAJ: [J{J + + S{S + l) - E^] , 



C.2. Spin-spin coupling Hamiltonian 



{ASj:njM\Hss\ASEnjM) = |a [se^ - s{s + 1)] , 



C.3. Rotational Hamiltonian 



{KST!q!JM\ H,ot\^SmjM) = 5(5s'Efe'n [J{J + 1) - + S{S + 1) - E^] 

-2BE,.«(-i)-— ( : ^ ) ( : I ) 

X y J( J + 1)(2J + 1)5(5 + 1)(2S + 1). 



C.4. Centrifugal distortion Hamiltonian of order 4 



( AS^'n'JM\H^^^ |A5EJ^ JM) = -D^s'S'^j^'n [^(^ + l)-n'^ + S{S + 1) - E^] 



J{J + 1)(2J + 1)5(5 + 1)(25 + 1) - 2L> Eg=±i(-l)''"^'+^~''' 

VJ(J + 1)(2J + 1)5(5 + 1)(25 + 1) 



J ^ J \ ( s is' 
-n' q n j y -T.' q e 

X [2 J( J + 1) - 02 _ (o')2 + 25(5 + 1) - E2 - (E')2 



C.5. Centrifugal distortion Hamiltonian of order 6 



AS^'n' JM\Hf^\AS^njM) = i?(5E'Efe'n [JiJ +1)-^'^ + 5(5 + 1) - E^]^ 



-31 - 



J1J\ S1S\ 

.q=±i ^w'l." y-n q n" J y-j: q s" ) 

xJ{J + 1)(2J + 1)S{S + 1){2S + 1) [3J(J + 1) - 2^2 - (J2")2 + 35(5 + 1) - - (S")^] 

' J 1 j \ f S 1 s\ 



-n' q n 



; + 1)(2J + l)S{S + 1)(25 + 1) 



[J( J + 1) - + 5(5 + 1) _ s2 



+ [ J( J + 1) - {n'f + 5(5 + 1) - (E')'] ' 

+ [J( J + 1) - + 5(5 + 1) _ S2] [J( J + 1) _ (J^/)2 + 5(5 + 1) _ (s/)2] 



+4 



J 1 J 



S 1 5 
-S' q S 



J(J + 1)(2J + 1)5(5 + 1)(25 + 1) 



(C5) 



C.6. Spin- rotation interaction Hamiltonian 



( ASJ:'n'JM\H,,\ASY.njM) = 7^s's<Jn'n [^^ - S{S + 1)] + 7Eq=±i(-l)^""'+^"^' 
/ J 1 J \ f 5 15' 



-n' q n 



-S' q E 



VJ(J + 1)(2J + 1)5(5 + 1)(25 + 1) 



(C6) 



C.7. Zeeman Hamiltonian 



(A5S'J'Mjfi'|i7z|A5E JMjfi) = /xbSq E<,=o ^^(-l)^-^'-^-"' [(2J' + 1)(2J + 1)]^/^ 



J' 1 J 



J' 1 J 



-M M M -o' g n 
s 15' 



5(lA5se' + (55 + 5'r)(-l 



,5-E' 



-E' g E 



[5(5 + l)(25 + l)] 



1/2 



griJ'BBoM6jjf6j:T.'6a,n' 



(C7) 
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Fig. 11. — Vector diagram for Hund's case (a) (left panel) and Hund's case (b) (right panel). 
Herzberg (1950). 
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